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ABSTRACT 

In this paper we develop an analytical theory for the interaction between waves and currents 
induced by breaking waves on time-scales longer than the individual waves. We employed the wave- 
averaging procedure that is commonly used in tlie literature. The near-shore zone is often 
characterized by the presence of breaking waves. Hence we develop equations to be used outside the 
surf zone, based on small-amplitude wave theory, and another set of equations to be used inside the 
surf zone, based on an empirical representation of breaking waves. Suitable matching conditions are 
applied at the boundary between the offshore shoaling zone and the near-shore surf zone. Essentially 
we derive a model for the interaction between waves and currents. Both sets of equation are obtained 
by averaging the basic equations over the wave phase. Thus the analytical solution constructed is a 
free vortex defined in both shoaling and surf zones. The surf zone solution is perturbed by a 
longshore component of the current. Thus the presence of the rip current cell combined with the 
longshore modulation in the wave forcing can drive longshore currents along the beach. The 
outcome, for our set of typical beach profile, is a description of rip currents. 

KEYWORDS: Wave-current interactions, surf zone, shoaling zone, matching conditions, wave- 
averaging, rip currents, radiation stress. 



I. INTRODUCTION 

Wave transformation in the surf zone is the dominant factor in the hydrodynamics of the nearshore 
circulation, and consequent sediment transport. A global description in terms of the spatial variation of such 
quantities such as wave action, wave action flux and wave radiation stress are the driving entities we have used 
to describe the generation by waves of mean currents in the nearshore zone. Studies on the interaction between 
waves and currents span nearly half a century. Mainly driven by a combination of engineering, shipping and 
coastal interests, there has been much research on shoaling nonlinear waves, on how currents affect waves and 
how waves can drive currents. This last aspect is the main concern in this paper. 

The basis for this subject was laid down by Longuet-Higgins & Stewart(1960, 1961), who analyzed the 
nonlinear interaction between short waves and long waves (or currents), and showed that variations in the 
energy of the short waves correspond to work done by the long waves against the radiation stress of the short 
waves. In the shoaling zone this radiation stress leads to what are now known as wave setup and wave setdown, 
surf beats, the generation of smaller waves by longer waves, and the steepening of waves on adverse currents, 
see Longuet-Higgins & Stewart (1962, 1964). The divergence of the radiation stress was shown to generate an 
alongshore current by obliquely incident waves on a beach (Longuet-Higgins 1970). 

The action of shoaling waves, and wave breaking in the surf zone, in generating a wave-generated 
mean sea-level is well-known and has been extensively studied, see for instance the monographs of Mei (1983) 
and Svendsen (2006). The simplest model is obtained by averaging the oscillatory wave field over the wave 
phase to obtain a set of equations describing the evolution of the mean fields in the shoaling zone based on 
small-amplitude wave theory and then combining these with averaged mass and momentum equations in the 
surf zone, where empirical formulae are used for the breaking waves. These lead to a prediction of steady set- 
down in the shoaling zone, and a set-up in the surf zone. This agrees quite well with experiments and 
observations, see Bowen et al (1968) for instance. However, these models assume that the sea bottom is rigid, 
and ignore the possible effects of sand transport by the wave currents, and the wave -generated mean currents. 
Hydrodynamic flow regimes where the mean currents essentially form one or more circulation cells are known 
as rip currents. These form due to forcing by longshore variability in the incident wave field, or the effect of 
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longshore variability in the bottom topography (Kennedy 2003, 2005 ,Yu & Slinn 2003, Yu 2006 and others). 
They are often associated with significant bottom sediment transport, and are dangerous features on many surf 
beaches (Lascody 1998 & Kennedy 2005). 

There is a vast literature on rip current due to wave-current interactions, see the recent works by 
(Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu 2006, Falques, 
Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others) and the references therein. 
Our purpose in this paper is to extend the wave-current interaction model examined in [5] to the quadratic depth. 

In section 2 we record the usual wave-averaged mean field equations that are commonly used in the 
literature. In section 3 we introduce a description of the rip current formation and examine the consequences for 
both shoaling and surf zones. Then in section 4 we employ section 3 to extend our previous results see 
Osaisai, E. F (2013), now to include the rip currents behavior in a quadratic depth profile. We conclude with a 
discussion in section 5. 

II. MATHEMATICAL FORMULATION 

2.1 Wave field 

In this section we recall the wave-averaged mean flow and wave action equations that are commonly 
used to describe the near-shore circulation (see Mei 1983 or Svendsen 2006 for instance). We suppose that the 
depth and the mean flow are slowly-varying compared to the waves. Then we define a wave -phase averaging 

operator < / >= / , so that we can express all quantities as a mean field and a wave perturbation, denoted by a 
"tilde" overbar. For instance, 

C = C+C- (1) 
where £ is the free surface elevation above the undisturbed depth h = h ( x ) . Then outside the surf zone, the 
representation for slowly- varying, small-amplitude waves is, in standard notation, 

C, \x,t) : a cos 0. (2) 
Here a = a(x.t) is the wave amplitude and = 0{x,t) is the phase, such that the wavenumber k , frequency 
fi are given by 

k = (k,i) = v $ a = -0 r (3) 

The local dispersion relation is 

Q = co + k.U , co 2 = gK tanh kH (4) 
where k 2 = k 2 + I 2 . 

Here U (x,t) is the slowly- varying depth-averaged mean current (see below), and H (x,t) = h(x) + £ (x,t) . 
To leading order, the horizontal and vertical components of the wave velocity field are respectively 

k cosh k (z + h ) _ sinh k ( z + h) 

u : — a Q cos , w : a Q. sin . (5) 

k sinh kH sinh kH 

Importantly, note that we have ignored here any reflected wave field, which is assumed to be very 
weak when the bottom topography is slowly varying. The basic equations governing the wave field is then the 
kinematic equation for conservation of waves 

k f + V co = 0, (6) 

which is obtained from (3) by cross-differentiation, the local dispersion relation (4), and the wave action 
equation for the wave amplitude 

A, + V -(c s A) = 0. (7) 

Here A = Elco , where E = ga 2 72 is the wave energy per unit mass, and 
c = Vj - co = U + c kl k ,(c = dcoldK) is the group velocity. Using the dispersion relation (4) in (6) we 
get 
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k t + c g • V k = -V „Q, (8) 
Q 1 + c i 'VQ = Q J . (9) 
Here the subscript "e" denotes the explicit derivative of Q (k , x , t ) with respect to either x or t , when the 
wavenumber k is held fixed. In this water case these explicit derivatives arise through the dependence of fi on 
the mean height H and the mean current U . 

2.2 Mean fields 

The equations governing the mean fields are obtained by averaging the depth -integrated Euler 
equations over the wave phase. Thus the averaged equation for conservation of mass is 

H t + V ■ (HU ) = 0. (10) 

Here H = h + £ where h = h(x) is the time-independent undisturbed depth. For the velocity field we 
proceed in a slightly different way, that is we define 

u = U + u , (11) 
where we define U so that the mean momentum density is given by 



M = HU 



= < f udz >, (12) 



But now we need to note that u does not necessarily have zero mean, and that U and u are not necessarily 
the same. Indeed, from (11) and (12) we get that 



u = U + < u > , and < f u 

J-i, 



dz >= 0. 



But u = u + O (a 2 ) , so that < u > is O (a 2 ) and it follows that, correct to second order in wave amplitude, 

E k 

M = Hu + M , where M = -H < u >=< Q u{x,0, t) >= . (13) 

C K 

The term M w in (13) is called the wave momentum. 

Next, averaging the depth-integrated horizontal momentum equation yields (Mei 1983) 

{HU ) t + V .(HUU )=- V. < f u u +pIdz> + <p(z=-h)>Vh. 

■ -h 

Next an estimate of the bottom pressure term is made by averaging the vertical momentum equation to get 

< p(z = -h) > - g(£ + h)= V .< [ wudz > + < \ w dz > . 
For slowly-varying small-amplitude waves, the integral terms on the right-hand side may be neglected, and so 
< p(z = -h) > » g (<£" +/?). Using this in the averaged horizontal momentum equation, and replacing the 
pressure p with the dynamic pressure q = p + (z - £ ) yields 

HU ( + V .( HUU ) = -V.S-gHV£ (14) 

( g ~ i 

where S =< \ [uu + ql ]dz > - < — f > 7. (15) 

l-i, 2 

Here 5 is the radiation sttess tensor. In the absence of any background current, so that U is O {a 2 ) , 
we may use the linearized expressions (2, 5) to find that 

S *c g k— +E[—--]I. (16) 
co c 2 

where the phase speed c = ml ' k , correct to second order in the wave amplitude. 
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2.3 Shoaling zone 

These equations hold in the shoaling zone outside the surf zone (defined below). In summary, the 
equations to be solved are that for the conservation of waves (6) combined with the dispersion relation (4), the 
wave action equation (7), the averaged equation for conservation of mass (10) and the averaged equation for 
conservation of horizontal momentum (14), where the radiation stress tensor is given by (16). In this shoaling 
zone, we assume that wave amplitudes are small, and that there is no background current. Then all mean 

quantities are O ( a 2 ), and in particular we can systematically replace H with h throughout these equations. 

Next we shall suppose that h = h(x) depends only on the offshore co-ordinate x > , where the 
undisturbed shoreline is at x = defined by h (0) = . Further, in the near-shore region, including all of the 
surf zone, we shall assume that h r > . Then we seek steady solutions of the equation set in which all variables 
are independent of the time / , and are also independent of the transverse variable >• . It then follows from the 
mean mass equation (10) that HU is constant, and since H = at the shoreline, it follows that we can set 
U = everywhere. Then in the dispersion relation (4) CI = co . From the equation for conservation of waves 
(6) we see that the frequency co and the transverse wavenumber / are constants, and the the offshore 
wavenumber k is then determined from the dispersion relation (4). As is well-known, it then follows that as 
H —> , I k I— > oo , that is the waves refract towards the onshore direction, where we assume that the waves 

are propagating towards the shoreline so that k < . The wave action equation (7) reduces to Ec e is constant. 
Near the shore, we can assume that the shallow water approximation holds and then c g w - ( gh ) " 2 , so that 

a 2 h " 2 w a~h l ' 2 , (17) 
where a is the wave amplitude at a location offshore where h = h.. The surf zone x < x, , h < h, = h(x, ) 
can then be defined by the criterion that h b is that depth where alh = A , defining an empirical breaking 
condition. A suitable value is A = 0.44 , see Mei (1983) or Svendsen (2006). 

The last step is to find the wave set-up C, from the mean momentum equation (14), which here 

becomes 

S x + gH C , = 0„ (18) 



1 



-)E. 



where S = — E C os 

c c 2 

Here <j> is the angle between the wave direction and the onshore direction, and S is the " xx " component of 
the tensor S . As h — > 0, <j> ~ Q, c g & c, S ~ 3E/2 , and we recover the well-known result of a wave set- 
down in the shoaling zone 

2 

CI 

Ah 



£ = 



a o h o 



Ah' 



(19) 



Here we have assumed that C, is zero far offshore. Note that the first expression for f does not need the use 
of the shallow water approximation, as shown by Longuet-Higgins and Stewart (1962). 

2.4 Surf zone 

In the surf zone < x < x b ,0 < h < h b , we make the usual assumption (see Mei (1983) for instance) 
that the breaking wave height 2 a is proportional to the total depth H , so that 

2 ZJ 2 



2a = yH 



E = 



(20) 



Here the constant y is determined empirically, and a typical value is y = 0.88. . To determine the mean 
height H = h + C, , we again use the mean momentum equation (14), but now assume that 



Issn 2250-3005 



|October| I2013 



Page 50 



Modelling the rip current flow- 



S = 3£/2 = r gH 2 /2 where r = 3y 2 /8 , so that 

h-h, 

THH + H (H - h) =0, so that H=H / + — , (21) 

(1 + D 

where the constant H = h b + £ b is determined by requiring continuity of the total mean height at x = x . 
Note that using (19) 

TT 1 2 I 1/2 I A 1 312 

H b = h b ~ a o h o /4h „ > 

and since H must be positive, there is a restriction on either the deep-water wave amplitude a or on the 
breaker depth h b , 

Note that the expression (21) is valid for any depth h ( x) , although in the literature it is often derived only for a 
linear depth profile h = ax . 

We are now in a position to determine the displaced shoreline x = x , defined by the condition that 

H = . That is, if h s = h{x s ) then H = (h- h s )/l +T , where 

h s = -Th b -(1 + T)C b , (23) 
Note that to use the expression (23) it may be necessary to extend the definition h (x ) into x < . For instance 
for a linear beach, h = ax , this is straightforward, but for a quadratic beach profile, h = /?jc 2 , the extension 
for negative x should be h = -/fx 2 say. Note that from (19), 

2 , 1/2 

. _ a o h o 

* - . 3/2 ' 

4/2, 



and, on combining this with the condition (22), it follows that the shoreline recedes (advances), that is 
h < 0(> 0) when 



a o < r 

4^ 2 l + r^' 

r /;f 2 a 



, 5/2 ,, 5/2 , 0/1 , 

< \ /A„ • (24) 



b ' 

l + rv 4// - 

This anomalous result does seem only to have been recently noticed [5] [5] and [5]. Since there is an 
expectation that the shoreline should advance (see Dean and Dalrymple (2002) for instance), essentially it states 
that the present model is only valid for sufficiently small waves far offshore, defined by the first inequality in 
(24), which slightly refines the constraint (22). Next we can normalized the generalized bottom profile 

H =h-h\ + r. (25) 

so that it can be written in the form 

H 1 h h 

— = (—-—)• (26) 

h b i + r h b h b 

Thus all the depths profiles we have examined can be plotted as normalized functions. 
2.5 Linear depth 

Now for the linear depth, this is just, for our parameters 1/(1 + r )( xlx b - x Ix. ) where 

2 , 5/2 

x a. «„ 

— =-r + (i + r) ° ° (27) 

and so depends on the wave input allh^ and the ratio h Q lh b . The normalized plots are shown in 
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figure [1] below. It shows wave height for linear and quadratic beaches. Observe, in linear case the slope is, 
a/(l + T) , and hence is smaller than the undisturbed slope, a . The plot is a function of xlx. from x lx. to 

1 where we have evaluated xlx, . Hence the plots depend on these two parameters, and we choose say 



a /h = 0.1 and h Ih = 1/4 (gives x < ), which gives x lx. 



-0.2. 



2.6 Quadratic depth 

The general bottom depth profile is given by equation (25). Our normalized parameters are 
1/(1 + T)(x I x I lx 2 b - (xjx b ) 2 ). Thus equation (27) is the same as that above for the linear depth case, except 
that the left-hand side becomes 

5/2 



-r + (i + d 



a K 

4h 2 h 5 ' 1 



Here the depth is + /? x~ where the sign is for x < , or x > which can be written as p I x I x. It follows 
that for the same parameters the right-side is again - 0.2 , and so x lx. = -0.45. Figure [1] shows that they 

agree at x = as the expression (25) already shows, but that the linear depth gives a greater setup in x > , 
but is weaker in the region x < . 
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Figure 1 : Plot of normalized wave height given by equation (25) for linear and quadratic profiles which 
depends on the ratios a /h = 0.1 and h./h. = 1/4. The values of x ranges from -0.2 < xlx. < 1 for the 
linear depth and - 0.45 < xlx h < 1 for the quadratic depth profile. Thus the graphs are plotted in the range 

- 0.45 < xlx, < 1 . 



III. THE GENERAL DESCRIPTION FOR THE RIP CURRENT FORMATION 

Here we consider a steady-state model driven by an incident wave field which has an imposed 
longshore variability. The wave field satisfies equation (7) which in the present steady-state case reduces to 

88x(Ec g cos 0) + ddy(Ec g sin 0) = 0. (28) 
Here we again assume that h = h{x) and that consequently the frequency m and the longshore wavenumber 

/ are constants, while the onshore wavenumber K is then determined from equation (4). We have the wave 
energy E of the form 



E = E (x)co& Ky + F (jc)sin Ky + G (x), 
where the longshore period 2 nlK is imposed. These equations in the shoaling zone yields 
(E n c cos 0) + KF „c sin = 

(F„c cos 9) - KE „c sin = 
(G.c cos 0) = 0. 

g ' x 



(29) 

(30) 
(31) 
(32) 



Issn 2250-3005 



|October| I2013 



Page 52 



Modelling the rip current flow- 



on collecting terms in cos Ky , sin Ky and the constant term, which form three equations for E , F Q 

and G . Equation (32) easily yields that G c cos 6 = constant . In shallow water, we may approximate by 

putting c ~ gh 12 and cos » 1 , so that then G ~ constant Ih 12 . For the remaining equations we can use 

Snell's law, sin 61c = sin 6 b lc, = a. (the constant value, here evaluated at the breaker line), and the shallow- 
water approximation to get that 

{( E c) x )/c} x + K 2 a 2 b E c 2 = 0, (33) 
while although F satisfies the same equation, once E has been found, then F is given by either (30) or 
(31). In practice, Kc « 1 and so approximately we can assume that (E ,F Q )c ^constant , the usual 

shallow- water expressions. Note that here c « J gh . In the surf zone, the expressions E (x), F ( x), G (x) is 
determined empirically. 

Once the expression (29) has been determined, we may then substitute into the expressions (30,31 & 
32) to obtain the radiation stress fields. Our aim here then is to describe how steady-state rip currents are forced 
by this longshore modulation of the incident wave field, especially in the surf zone. 

The forced two-dimensional shallow water equations that we use here are characteristic of many 
nearshore studies (Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu 
2006, Falques, Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others). Then, 
omitting the overbars as before, then equations (14) in the present steady-state case reduce to 

H[UdUdx + VdUdy] = -g Hd^dx - [rj, (34) 
H [UdVdx + VBVdy] = -g Hd^dy - [r J, (35) 
where the stress terms are defined; 

r r = 8S n 8x + 8S 12 8y and r y = 8S 21 8x + 8S 22 8y . (36) 

Next we observe that equation (10) can be solved using a transport stream function y/ (x, y) , that is 
1 1 

U = - — 8y/8y and = — 8y/8x, (37) 
H H 

Next, eliminating the pressure, we get the mean vorticity equation 
Q. Q r 

yf , (— ) -w y (— ) ( = [— ] , " [— ] « (38) 
H H H H 

where Q is define as 

V V v 

Q = V x -U = (—), + (—) (39) 

H H 

We shall solve this equation (38) in the shoaling zone x > x, and in the surf zone x < x. , where as 

before x = x b is the fixed breaker line. It will turn out that the wave forcing occurs only in the surf zone, but 
continuity implies that the currents generated in the surf zone must be continued into the shoaling zone. 

3.1 Shoaling zone 

In x > x, we shall assume that H a h as £ is O {a 2 ) . Then we shall use the expressions [30 ,31] 
to evaluate the radiation stress tensor. For simplicity, we shall also use the shallow-water approximation that 
c « c « J gh , and so we get that 
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S n = £(cos 2 + -),S a = 5„ = E sin cos , S 22 = E ( sin 2 + -) (40) 
2 2 

These expressions are in principal known at this stage, and so we can proceed to evaluate the forcing 
term on the right-hand side of (38). To assist with this we recall Snell's law 

sin = -Jh -Jh^ sin 6 b 

where h. and are the water depth and incidence angle at the breaker-line. Now the energy equation (28) has 
the approximate form 

(Eccos 6) x + (Ecsin 0) = 0, 
and using Snell's law, this can be written as 

(E cos 2 0) t + (E sin cos 0) + E — = 0, 

c 



1 c x 
and so r = — E — E — -. 

2 ' c 

We can also deduce from (28) that 

( E sin 9 cos 9 ) x + {E s in 2 # ) „ = , 



and .vo r = — E . 

y 2 " 

We can now evaluate the right-hand side of (38), and find that its identically zero, 

T « r v 

[—]-[—], = 0. 

h h 

Thus in the shoaling zone there is no wave forcing in the mean vorticity equation, although of course there will 
be a mean pressure gradient. However, this does not concern us since here our aim is to find only the flow field. 
Note that the result that there is no wave forcing in the vorticity equation does not need the specific form (29), 
and is based solely on the steady-state wave energy equation (28). The specific form (29) is only used in the surf 
zone. 

With no forcing term, the vorticity equation (38) can be solved in the compact form, noting that we 
again approximate H with h , 

— =F(yy). (41) 
h 

But here F {y/ ) = from the boundary conditions in the deep water as x — > oo , where the flow field 
is zero. Thus our rip current model has zero vorticity in the shoaling zone. It follows that we must solve the 
equation 



1 1 

Q = {—yr ) + (-y, ) = 0, 
h h 



(42) 



in x > x, . Since h = h(x) we can seek solutions in the separated form 

yr = X (x)Y(y) 

with the outcome that 



X 

(— ), 
h 



K 2 X 



h 



0, Y + K Y = 0. 



(43) 



(44) 



We note the separation constant K =2 ttIE must not be zero, and is in fact chosen to be consistent with the 
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modulation wavenumber of the wave forcing. Without loss of generality, we can choose 

Y = sin Ky . 



(45) 



For each specific choice of h{x) we must then solve for X (x) in x > x, , with the boundary 
condition that X — > as x — > oo . We shall give details in the following subsections. Otherwise we complete 
the solution by solving the system (38) in the surf zone, and matching the solutions at the breakerline, x = x b 
where the streamfunction y/ must be continuous, and in order to have a continuous velocity field we must also 
have that y/ t is continuous. 

3.2 Surf zone 

To make sense of wave forcing, we assume that the expression (29) holds in this region. The functions 
E Q (x), F (x), G (x) are then determined empirically. To determine the wave forcing term in the mean 

vorticity equation (38) we shall assume that = « 1 so that, on using (??) and (40) we get that 



t = — E 
2 



1 

T = — E . 

2 * 



Then (38) now becomes, where we again approximate H with h(x) 



1 y/dxdQ. dy - 8 y/dydQ-dx = — — + 



Eh. (h E ) 



2h- 



(46) 



where here Q = Cl/h is the potential vorticity. Since the wave forcing is given by (29), that is 
E = E (x)cos Ky + F (x)sm Ky + G (x), 



(47) 



we observe that the unmodulated term G ( x ) plays no role here at all, although of course it will 

contribute to the wave setup. In order to match at x = x, with the expression (45) for the streamfunction in the 
shoaling zone, we should try for a solution of (46) of the form 

iff = F (x) sin Ky +G(x), in x < x b . (48) 



The matching conditions for the streamfunction and velocity field at the breakerline x = x require 



that 



F(x b )= X(x h ),F x (x b ) = X x (x h ),G(x = x b ) = 0,GJx h ) = 0. 
The expression (48) yields 

Q = F sin Ky + G 



(49) 



where F and G are differential operators where they are defined as; 

F K 2 F 

F = (— ) x 

h h 



G = Z , Z = 



h 



(50) 
(51) 



From equation (46) we get a set of three equations that are used to determine the rip-current flow field 
in the surf zone. These are namely; 



F F 
F —- F(— ) = 0, 
h h 



(52) 
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F G (h m F n ) 

(G—-F (—),.) = (53) 

h h h 

h m E ) 

( d^=0. (54) 

h m 

Equation (54) gives E a :l/ h" 2 , which is an unacceptable singularity as h —> . Hence we must infer 
that in the surf zone at least, E = . The first of the three equations, that is (3.2a) suggests that 

F 

— = CF whereCisac onstant , (55) 
h 

and the second (3.2b) yields that 

G (h m FJ 

F(CG -(—),) = - (56) 

h h 

The boundary conditions at x = where h = are that both mass transport fields t/ , V should 
vanish, that is from (37) i// = constant and ly Ih = , which implies that 

G 

F = F = , G = constant , — — = , at x = . (57) 

h 

As above there are also the matching conditions for both F and G separately at the breakeiiine, that is for F 
we have that 

= , at x = x 

F(x b ) X(x h ) 

where we note that here the right-hand side is a known quantity, depending only on K and x b . Next we see 
that equation (55) reduces to 

2 

F K F 

(—) = ChF (58) 

h ' h 

Together with the boundary conditions at x = 0, x = x. this is essentially an eigenvalue problem for 
F (x ) with eigenvalue C . In general it is solved approximately since we shall assume that Kx « 1 . Once 

F (x) is known we can solve (56), together with the appropriate boundary conditions to get G ( x ) to complete 
the solution. 

Note that the amplitude of F ( x ) is an arbitrary constant in this solution, and so we can fix it by 

specifying its value at x = x, say. Indeed the solution we have constructed is essentially a free vortex defined 

by X (x)sin Ky in the shoaling zone x > x , and F (x)sm Ky in the surf zone x < x , perturbed by a 

longshore component G ( x ) in the surf zone. Note that in the presence of the wave forcing, both F , G are 
non-zero, see (56). It is significant that unlike the longshore currents considered in the basic state which depend 
on an ad hoc frictional parametrization, the presence of the rip current cell combined with the longshore 
modulation in the wave forcing can drive a longshore current. 

IV. APPLICATION TO A QUADRATIC DEPTH PROFILE 

Osaisai, E.F (2013) examined the behavior of the rip-current for the case h = a x. We here extend the 

case for which h = p x ~ . In the basic state we show that the linear depth gives a greater setup in x > , but is 
weaker in the region x < , where both depths agreed at x = . 
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4.1 Shoaling zone 

Now we let h = fix 2 and then equation (44) now admits the differential equation of the form 

xX " (x) - 2 X ' ( x) - K 2 xX ( x) = 

whose solutions are explicitly the exponential functions of the form 

X (x) = C.e & (Kx - 1) + C,e~ & (Kx + 1). 



(59) 



From the behavior of the solutions as 



we see that C , = and so 



X (x) = C 2 e~ r (Kx + 1) . 
Again note that the constant fi does not appear in this solution. 

4.2 Surf zone 

Similarly, let ft = fix 2 in (58), so that proceeding as for Osaisai, E.F (2013) we obtain the equation 



(60) 



K F = A n Xx , 



(61) 



where now A = C fi 2 , and as x — >• , F » A Q x^ . The solutions of the homogeneous equations are 

known, these are j = e (Kx + 1) and j = e (Kx - 1) where /, and j 1 are linearly independent 
solutions. The general solution is given by 

F(x) = A(x)j, + B(x)j. + C.j. + C.L, (62) 



where the Wronskian of j\ and j 2 is W = 2 K 3 x 2 , and 

A(x) = - — 5 — f j 2 x 5 dx 
2K Jo 



12 3 



A A >- 5 A Ajc 

fi(x) = j r x dx ~ — . 

2K Jo 12 K 

Thus A and B vanish at x = as x 5 and so C , = C 2 . Also the normalization of F as jc — » gives 
C, = 3A () /2 ^ 2 . 

Now to apply the boundary condition at x = x. which yields A , we simplify the calculation by first 

approximating A and B as above for small x . The outcome is that 

, i 
, Ax 

F m A (x + ). 

36 

The boundary condition at x = x. again yields an explicit equation for A which can be simplified by the 
assumption that Kx « 1. Thus A scales as x^ 6 which may not be so good an approximation. In spite of 
that the RHS may be approximated by - K 2 x\ . Hence to leading order 



Ax. = -12 . 

b 

This leads to the simple expression 



(1 -— ), 



(63) 



(64) 



3x u 
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and evidently F » at x = x b . 




Figure 2: Plot of F (x)IA against xlx. , where A is arbitrary as given by equation (64). Observe there is no 
dependence on the slope /? but only a weak dependence on K . The value of F (x) also reaches a maximum 
value of F(x)/A « 0.7. This shows that irrespective of h(x), the maximum value F (x) can admit is 

a 0.7 . 



that 



As in [5] we can now add correction terms, letting s = 1+ Xx b IYl . Expanding A and B we find 



1 2 

F (x)= A (x J -— K"x" - K"x" + A[ x" + — K V ). 

4 72 36 9 



As before, we now find the leading order term for e 



e = 1 + Zx l\2 



Finally we get that 



x Kx 



F (x) = A x (1 K x- - 

4 3xt 3x1 



(65) 



Next, as before, we need to solve for G (x) from (56). As above, we approximate F = A {) .* 3 , and also we use 
the empirical expression F = y 2 h 2 /8 , see equation (20). Thus we get 



Z - — Z - C B x Z = 

XX X ' 

x 8A„ 



(66) 



Letting u = x 3 we get that 



9Z - AZ 



5gy 2 /3 2 



SA a u 3 

As before the dominant balance in the particular solution is between Z and the right-hand side, so that 
Z fa constant u "' 3 . We find that 



Z = 

p 



at 1 n 1 2 

45 gy P x 
16 A n 



(67) 
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As in the linear depth case, this have the form as X < 



J X u V X u 

Z , = C , sin + C , cos 



Here we recall that X = - X . The total solution is then 

Z = Z + Z, . 

The boundary condition at u = gives C 3 = 0. Again imposing the boundary conditions that Z = at the 
breakerline x = x gives, 



Ac- n 2 2 2 2 

45 gp r *t_ * 

16 A„ *, 2 



2x- 
sin 1-^=1). 



V— L I 

3 ] V 3 ^ 



(68) 



Finally we get G from G = hZ and G = at x = x 



45 p x x 
Cr = ( r + 



1 



2x 1 
cos [ — i ] 



16 A 5x 4 5 2V3~sin [2/ V3] 5 tan [2/ ^3] 



(69) 




Figure 3: Depicts the plots of normalized Z ( x ) and G ( x ) given by equations (68) and (69) where each is 

normalized by 45 g p 2 y 2 x~/16 A and 45 g p 2 y 2 x 5 b l\6 A respectively with A arbitrary. We observe as in 
[5], here too as depicted in the figure, there is a small region of reversed flow near the breaker line. 



The combined expressions (60, 64, 69) complete the solution, where we recall that the constant C is 
given by (63) (since X = C P 2 ), or their respective higher-order corrections. Now the amplitude of F (x) at 
x = x b is given by 

F(x b ) = C 2 e~ Kx (Kx +1). (70) 
On using the approximation Kx « 1 , and the approximate expression (64), this reduces to 

2A„jc, 3 
F(x h ) = C 2 . 
3 

The rip-current system contains a free parameter A or its equivalent. We choose to define this free parameter 
to be the value of F (x b ) and normalize the full solution by this value. Thus we get from (43, 45) in x > x. , 
and (48, 64, 69) in x < x. that the normalized streamfunction y/ n is given by 



_ 
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X (x) 

if/ = sin (Ky ), for x > x 



V n = 



F(x) 



G(x) 

sin Ky + R , for < x < X 



F(x h ) G(0) 
Here again R = G (0)/ F (x. ) is a free parameter. From (64, 69) we find that here 

135 gj3 } r 2 x; l l l 



32 A 2 2V3 sin [2/ tan [2/ 

Note that again i? < , and that I R I increases as the wave forcing increases, or as the curvature 
increases, or as the depth fix * at the breaker line increases. In order to estimate typical values for R we again 
note that from (64) the longshore velocity field in the "sin {Ky ) "-component scales as V = A //3x b , while 
the longshore component then scales with RV . Taking account of the actual numerical values in the 
expressions given above, we find that a suitable values are R « - 0.1 . Plots of y/ i: are shown in figure [4 & 5 
] for same values of R as in the linear case, and again Kx = 0.2 . 



(71) 
(72) 

(73) 




0.5 1.0 1.5 2.0 2.5 3.0 



1.0 1.5 2.0 2.5 3.0 



Figure 4: Plot of the rip current streamlines for a quadratic depth profile, given by equation (72) where F (x ) 
and G ( x ) are equations (64) and (69) respectively for if = - 0.02 in the left panel and R = - 0. 1 in the right 

panel. 




Figure 5: As for figure 4 but R = -0.5 in the left panel and R = -2 in the right panel. 
Overall these plots show the same kind of behaviour as those for the linear depth profile see [5]. 
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However, the major difference is that the flow in the surf zone is rather weaker, and so the vortex centre is 
slightly further offshore. 

V. CONCLUSION 

We described qualitative solutions for rip currents which are essentially free vortices in both zones. 
The free vortex in the surf zone is perturbed by a longshore modulation in the wave forcing. Rip current cell 
combining with the longshore modulation in the wave forcing can drive longshore currents along the beach. 
Thus the dynamics of the shoaling zone is only dependent on the state-state wave energy equation.The wave 
forcing in the surf zone sets the wave activities different from those of the shoaling zone. To determine wave 
forcing in the mean vorticity equation we assume that the wave angle becomes smaller. We also note here that 
the component of the radiation stress in the y momentum remains unchanged across the entire flow domain. 

This shows that it is only the x component of the radiation stress that play a leading role in the wave forcing. 
However, wave forcing encountered in the surf zone has an unmodulated term that does not play a role in the 
vorticity equation but only contribute to wave setup. To ensure continuity of the streamfunctions in the shoaling 
zone we match the solution at the breakerline by a matching condition with appropriate boundary conditions. 
Thus the rip currents solution in the surf zone is provided by the terms in the matching condition. The terms in 
the matching condition has a cross-shore width and a modulated longshore component. The cross-shore width 
was determined by the application of perturbation method and variation of parameter. It would to interesting to 
examine the effect of friction on the rip currents. 
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